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The interaction of a three-level A-system with a Markovian environment is analyzed. A coherence 
vector formalism is used and a vector form of the Lindblad equation is derived. Depolarization, 
dephasing, and different amplitude damping channels acting on the subspace of the ground states 
\ of the A-system are investigated. Analytic expressions of the concomitant optical susceptibility 

■ are given in the limit of a small probe field and a strong control field, and we analyze how the 

' different channels influence the optical response. For a given group velocity, the absorption depends 

^S) I non-trivially on the type of open system interaction and even gain can be obtained. 
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INTRODUCTION 



A medium that is opaque to radiation on one of the transitions of a A-system can be made transparent by addressing 
a laser field on the second transition, a phenomenon known as electromagnetically induced transparency (EIT) [H, Q- 
Associated with this transparency is an extremely small and dynamically controllable group velocity j3||. The pulse 
. can be brought to a complete stop and stored as an atomic coherence. Most impressive with respect to storage time, 
^ ' is perhaps the storage of a light pulse in a solid for a period greater than one second The process can be reversed, 
O^. and the light pulse released with the original pulse information intact 0], making the EIT a strong candidate for a 
quantum memory Beside applications in quantum information processing, the narrow spectral feature of EIT is 
[ an important tool in spectroscopy Q, high precision magnetic sensors Q, and atomic clocks [9l.[lo|. 
J> ■ The dark state of the A-system is responsible for the above coherence properties Open system effects suppress 
CNj ' the dark state of the system and thus sets the ultimate limit on the storage time of light, the residual absorption, and 
•/^ ' the spectral resolution of the transparency. There are tricks that one can play in order to increase the efficiency of this 
coherence effect. In, for example, Ref. Q the effective decoherence rate was decreased by a method called quantum 
bang-bang dynamic control, inspired by refocusing techniques in nuclear magnetic resonance fl^. Furthermore, there 
OO is a trade-off between the delay of a pulse and the bandwidth In order to increase the bandwidth, spatially 

dispersed Hght beams in non-isotropic media can be used fl3]. 

An electromagnetic field interacting with matter gives rise to a matter polarization. The Hnear part of the polar- 
ization field P — P{E) can be written in terms of the electrical susceptibility x — P/iEeo), with eq the electrical 
permittivity in free space. The real and imaginary parts of the susceptibility, describe the phase velocity of the 
electrical field and the absorption or gain of the field, respectively. The macroscopic polarization field of an isotropic 
medium consisting of N atoms per unit volume is 



P= (P> =A^(p> =A^Tr[pA], (1) 

where p is the atomic density operator and /i = e f is the dipole operator with electric charge e and position f . 
Therefore, if p is known as a function of the optical field, the optical response is too, as it is given by the coherences 
of p. 

In this paper, we investigate various types of interactions between the A-system and a Markovian environment and 
how they infiuence the optical response, i.e., the susceptibility of the medium, via calculation of the density operator 
p in the asyniptotic time limit. We use the Gell-Mann matrices to reformulate the Lindblad- Kossakowski master 
equation [H, employing a coherence vector formalism . A convenient vector form of the master equation is 
found, in which the dynamics is rephrased as a vector equation using SU(3) vector products. The coherence matrix 
representation of the vector master equation is given as an inhomogeneous linear matrix differential equation, to which 
the considered asymptotic state solutions are independent of the initial state. 
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The linear optical response of a weak probe field interacting with the A-system can be determined analytically and 
an explicit form of the susceptibilities due to the dephasing, depolarization, and three different amplitude damping 
channels acting on the ground states of the A-system are given. Effects due to open system interactions are expected 
to be relevant in many experimental instances. 

The outline of the paper is as follows. Section II contains a general theory for a Markovian dynamics in a three- level 
system using coherence vector formalism. In Sec. Ill, we apply this framework to our specific A-system, undergoing 
an open system dynamics corresponding to the dephasing, depolarization, and amplitude damping channels. We are 
interested in the asymptotic solutions of the system, i.e., when t —>■ oo. From these it is possible to deduce the 
susceptibility of the system and open system effects are explicit. The paper ends with the conclusions. 



II. OPEN THREE-LEVEL SYSTEM 



We commence by studying a general three-level quantum system and the effects due to the interaction with an 
environment. For a Markovian environment, the density operator p may be assumed to satisfy the master equation 

dtp=^[H,p\+Cip), (2) 

where H is the Hamiltonian operator. The Liouville superoperator £ is linear, describes the coupling with the 
environment, and can be written on Lindblad form 

^(/5) = J2^^^P^l ~ h^'l^^P ~ 5^7k7k)' (3) 

k 

where 7k are Lindblad operators, and the subscript k denotes the open system channels. 

The coherence vector representation consists of a real-valued vector x (the 'coherence vector') and the generators 
of the appropriate algebra For a three-level system, we express the dynamics with respect to some Hilbert space 
basis in terms of the SU(3) matrix generators, i.e., the traceless Gell-Mann matrices 



Ai = 1 , A2 = 



A, = -1 , A4 - 





As = , As = 
\i J 

/O \ 

At = I -I 1 , As = ^ I 1 I . (4) 
These satisfy the fundamental relations 

Tr[ArAs] = 26rs, 
4 

{Ar,As} = -Srs + '2drst>^t, (5) 

where / and d are real- valued structure constants of the SU(3) algebra, being totally antisymmetric and totally 
symmetric, respectively, under permutation of the indices r, s, and t [l^. The values of the structure constants are 
tabulated in, e.g., Ref. [131 ■ The number 1/3 in Eq. ^ refiects the dimensionality of the quantum system, for a 
general N level system this should be replaced with The density matrix p and Hamiltonian matrix H can be 

written in terms of the A's as 

p = ^1+ X- \, 

H = hwol + huj-\. (6) 
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Here, 1 is the 3x3 identity matrix, we have introduced eight-component vectors ui and A associated with the 
Hamiltonian and Gell-Mann matrices, respectively, and we have used the standard scalar product ' • '. The length of 
the coherence vector a; is < |xp < 1/3, with \x\'^ = 1/3 for pure states and \x\'^ — for completely mixed states. 
By inserting Eq. ([6]) into Eq. |[2]) we obtain 



X \ ^ -[u}-X,x-\] + ]-Cn)+C{x-\), 
I 6 



(7) 



i.e., the dynamics reformulated in terms of Gell-Mann matrices. Any physical process giving rise to a change in 
the length of the coherence vector can be split into a Hamiltonian contribution and a traceless Lindblad operator 
contribution to the Liouville equation. Thus, the matrix form of the Lindblad operators 7k is written as 7k = 
gol + Sk ' where the symbol '=' stands for 'represented by'. In this paper we will not consider the contribution 
from 5o, but focus on the 'Lindblad vectors' [2ll |. 

Our goal is to put Eq. ([7]) on A-independent form. To this end, we introduce the following vector products 

" A /3 = frsterCtsPt, 

Ct*(3 = drstGrasPt (8) 

where ol — aqSq and f3 — jSqEg are arbitrary complex-valued vectors, {ei}^^i is a real orthonormal basis of R®, and 
we have used Einstein's summation convention (repeated indices are summed). The wedge product 'A' is a higher- 
dimensional vector product analog to the cross-product in three-dimensional space. The star product ' is a vector 
product that has no counterpart in three dimensions. 

Equipped with these tools we may now rewrite each term in the right-hand side of Eq. First, the Hamiltonian 
contribution reads 

-[uj-\,x-\]^2iujAx)-X. (9) 



Secondly, the middle term of Eq. ^ is readily evaluated as 

£(l)-^2z (SkAflD-A, 



(10) 



where each term g^ A is purely imaginary, and hence £{t) is real as required. Note that, for Hermitian 7k the 
corresponding vector g^k is real and thus g^^ A g^. — 0. Finally, by using the symmetries of the structure constants the 
last term of Eq. ^ takes the form 

Cix-X) = ^ { [fif^ A (s-k A a;) + g-k A (g-k A a;) 

k 

+ ^[(9k Aa;)*g^+gk*(a;Agi^) -2(g;jAgk)*a; + 3(gk*a;) Ag;j + 3gkA (a;*g;:)]} • A. (11) 

By identifying terms, we arrive at the full Lindblad master equation in vector form 

'2i 



2uj A 

i 

+ 4 



-9k A gk - 



9k A (gk A a;) + gk A (gk A x) 



k 

(Sk A a;) * g;j + gk * (a; A gk) - 2(g;^ A gk) * a; + 3(gk ★ a;) A gk + 3gk A (a; ★ gjj) | . 



(12) 



This equation describes the dynamics of a quantum system that interacts with a Markovian environment. It should 
be emphasized that since the Gell-Mann matrices have been eliminated, this expression is valid for any iV-dimensional 
quantum system provided the SU(N) structure constants are being used in the definition of the vector products A 
and 

Let us now proceed by representing Eq. I|12p in a coherence matrix form. We write x = XiSi, u) — LUiBi, and 
A = A^Cj, yielding 

M("'a; : a4°^ = 2frstuJs, 



2ijj Ax 

9k A (gk A a;) + gk A (gk A x) 
(gk A a;) * gj^ + gk * (a; A gk) - 2(g^ A gk) * a; + 
-F3(gk ^x)Agl + 2,g^ A (a; ★ g*) 
^9k A gk 



f>k: 



5k,r 



3/, 

sv ' 

(13) 
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FIG. 1: The A-system. The field on the a-b transition is a weak probe field ilb, and the field on the a-c transition the control 
field D,c. The upper level is broadened by the finite lifetime I/7. 

with 

'^k% = i9k,s9Kv ± ffk.sffk,!,), (14) 

being symmetric (+) and antisymmetric (— ) in the indices s and v. Furthermore, while C^^J^ is real-valued, C^J^^ is 
purely imaginary and thus vanishes for real- valued gj^. The antisymmetric matrix M^^^ corresponds to the Hamiltonian 
of the system. The second row of Eq. ifTSl defines the real and symmetric matrix g\^\ The matrix G^^'' and also 

the vector 6k vanish for real-valued g^, as they both are proportional to the matrix C'^tK The matrix g\^J^ is also 
real-valued but has no obvious symmetry in the indices r and t. 

The complete dynamics of the system, including open system effects, can be written as an inhomogeneous matrix 
differential equation [22| 



with the evolution matrix 



and vector 



X = Mx + 6, (15) 
M = m(°)+^(4+)+4-)) (16) 

k 



b = ^6k. (17) 

k 

For time-independent and invertible M the general solution of Eq. (flSl) takes the form [TtI. fT3| 

X = e*^*a;o - M^^b, (18) 

where xq — x{0) + M^^b, x{0) being the initial coherence vector. Note that the evaluation of the exponential of M 
is non-trivial, since M is not diagonaHzable in general. 



III. OPEN A-SYSTEM 
A. Dynamical quantities 

The A-system consists of two ground state levels coupled to an excited state via electromagnetic fields, as shown 
in Fig. [TJ The corresponding Hamiltonian is given by H{t) = Hq — fi ■ E(i), where Hq is the Hamiltonian of the bare 
atom, p, is the electric dipole operator, and E(i) is an oscillating electric field. Furthermore, let Si be the energy 
eigenvalues of Hq corresponding to the bare eigenstates \i), i = a, b, c. In a frame rotating with the external field, the 
Hamiltonian operator reads 



H = hSb\b){b\ + hSc\c){c\ - hnb\a){b\ -hnc\a){c\ + h.c, 



(19) 
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where we have used the rotating wave approximation and we have defined the complex Rabi frequencies in terms of 
the electric field E in the rotating frame as fli = (a|AN) ■E/(2?i) and detunings Si = uji — {£a —£i)/h, where uoi are the 
frequencies of the field. Defining the Gell-Mann matrices in terms of the ordered basis |c), |a)}, the parameters 
of the Hamiltonian matrix in Eq. |[6]) are 



Wo 



Tr[i/] 



5e^ — |ilf,| coS(/)fce4 — sin(/)f,e5 



sm <pce^ 



(20) 



where Qi — \ \ exp(i^i) and we have defined the two-photon detuning 6 = {Sb ~ Sc)/2 as well as the mean detuning 
A= {6b + Sc)/'2. The factor ^ in the definition of S and Ifli] has been introduced for notational convenience. 

We now specify the Lindblad operators that describe the different forms of open system dynamics. First, sponta- 
neous emission, or the decay channel, is always present in the A-system. It corresponds to the amplitude damping 
of the excited energy state \a) to either of the two ground states |6) or |c) with rates jb and 7c, respectively. The 
corresponding Lindblad operators and their matrix representation read 



lb 



~b\b){a\ 



7fc(A4 + 1X5), 



Ic = \/>|c> (a| = - V7^(A6 + ih), 



(21) 



and they define the natural fine width 7 = 7b -f 7c indicated in Fig. [TJ The subscripts b and c indicate the final 
quantum state of the atom after interaction with the environment through the decay channel. The Lindblad vectors 
are 



9b 
9c 



2 V^(*^4 - 

= |^/7^(e6 



ies), 
iej). 



(22) 



These vectors are complex-valued and therefore the matrix C^^^ in Eq. (fT4|) is non-vanishing and, accordingly, the 
decay channel contributes to the dynamics, Eq. (flSl) . with a non-zero b vector. 
The dark state of the A-system is 



\d) 



1 



,{n:\b)-ni\c)). 



(23) 



This is an eigenstate of the Hamiltonian, it is also unaffected by the decay channels defined by Eq. l(2T|) . The 
additional open system effects that we are primarily interested in are the ones that act on the two-dimensional subspace 
Vbc = {\b),\c)}, since these could potentially destroy the dark state and its associated phenomena. Furthermore, we 
consider the limit of a weak probe field and a strong control field, i.e., \flb\ ^ \^c\, hence the dark state will be close 
to |6). 

For a general two-level system a set of generic open system channels is given by the depolarization, dephasing, and 
amplitude damping channels [l^l- Here, we examine these channels acting on Vbc- Depolarization and dephasing are 
particular combinations of the Hermitian Lindblad operators 



% = Vv^{\b){c\ + \c)(b\) = ^Xi, 
Iv = \Afe(-«l^)(c| +«|c)(6|) = V^Az, 
% = V^(|&)(6|-|c)(c|) = v¥A3, 



(24) 



where the symbol iji is used to denote rates of the various channels acting on Vbc- The corresponding Lindblad vectors 
are real and given by 



The depolarization channel is isotropic, i.e. 
hand, is anisotropic and corresponds to rj^ 
b vector, since the matrix C'-"-' vanishes. 



9x 
9y 
9z 



'7x61, 
^62, 

vies. 



(25) 



corresponds to rj^ = r]y = r]^ — ri/3. The dephasing channel, on the other 
= rjy ^ 0,7]^ ^ 0. Neither depolarization nor dephasing contribute to the 
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In the limit of a small probe field, the dark state is close to \b), and we may expect a certain asymmetry between 
amplitude damping in the two directions b ^ c and c^h. To take this into account, we choose to study combinations 
of the two amplitude damping channels 



7c6 = V^I^X^I = ^\/^(Ai - «A2), 
where 7hc fiips |c) to |6) at a rate rycfc and vice verse for 7c6. The corresponding Lindblad vectors read 

96c = 5\/%^(ei +«e2), 



9c6 



/^(ei - 162), 



(26) 



(27) 



which are complex-valued and may therefore contribute to the b vector. 

The evolution matrix M and the vector b contain contributions of three different types: 

(i) Decay of the excited state is considered to be part of the basic setting of the system, viz., the decay channel 
is always present in the analysis from this point and onwards. The corresponding M matrix has the form 



^ + Gc^^ + Gc '' and the bdecay vector is non-zero since the Lindblad vectors of the decay 



-(+) 



channel are complex. As we will see below, 6d 



bfc -I- be points in a fixed direction in W' 



(ii) Dephasing and depolarization in Vbc corresponding to real- valued Lindblad vectors g^, for which 6k vanishes 
identically. In this case the evolution matrix is given by = G^."^' . 

(iii) Amplitude damping in Vbc corresponding to complex- valued g^. In this case, one has to evaluate G'\^\ G^ \ 
and bk. 

The complete dynamics of the system is described by the evolution matrix M and vector b, where we are sum- 
ming over all channels, as in Eqs. I|16|17p . These quantities take a simple form by applying the following similarity 
transformation 



T 



T T 
^1^1 + 62^2 



'6^5 



The transformed evolution matrix M' is given by 

M' = TMT-^ 

the submatrices of which are given by 



A = 



A C 
-GT B 



(28) 



(29) 



B 



C = 

sm <Pb 

Here, we have introduced the parameters 



V- 
V+ 
Tb 
Lc 

r 




— |r2c| cos( 
\Qc \ sin ( 
-\Qb\ cos cpb 
-\/3|rib| cos (j)b 



Qb\ sin <j)b 

flbl cos (pb 

Sic I sin (jjc 
V3\^c \ sin ipc 



-|Slh| cos 0h 
-\^b\ sin (j)b 
|ric| cos <j)c 



COS 



(30) 



(f?6c - ?7c6)/2, 

Vbc + Vcb + 4?7, 
iSb + ^+ [iri + rjcb + %)/2, 
i5c + l + {irj + rjbc + ri^)/2, 
2i6 + (r/bc + Vcb)/2 + 4t] + 2rj^ 



(31) 



7 



The Hamiltonian contribution M^^^' in the antisymmetric part of M' and corresponds to an infinitesimal rotation 
of the coherence vector in R*. The diagonal elements of M' are real and negative, and contain the open system 
rates. The only channel rate that is located off the diagonal of the evolution matrix is rj_ / VS, in submatrix A, which 
originates from the two amplitude damping rates in Vbc having different channel rates. 
By using Eqs. I|13p and (|28l) . the transformed vector 

b' = Tb, (32) 

takes the form 

b' = + -^64. (33) 

The first contribution to b' stems from the two amplitude damping channels and vanishes only if the two rates yy^c 
and -qcb are equal. The second term on the right-hand side of Eq. (|33|) corresponds to the decay of the excited state 
\a). Note that the dynamics depends on the sum of "fb and 7c, not of the relative rates. 



B. Jordan form analysis 

The evaluation of the exponential of M is non-trivial, since M is not diagonalizable in general. To deal with this 
fact one may resort to the Jordan normal form of M, i.e., M = SJS~^, where J has non-zero entries only on the 
diagonal and the super-diagonal, and S is an invertible matrix. Expanding the exponent in Eq. (flSl) . we can write 
the general solution as 

X = Se^'S-^xn - M-^b. (34) 

The spectrum of M is {sijf^i, Sj € C, where the real part of the eigenvalues originates from C and must be zero or 
negative in order to preserve positivity of the density operator. 

We now look at the eigenvectors of M. For our A-system these will in general be non-orthogonal. Higher-dimensional 
Jordan Blocks occur provided M has linearly dependent eigenvectors corresponding to degenerate eigenvalues. How- 
ever, it is a rare situation to have two or more linearly dependent eigenvectors and in all the cases we investigate in 
this paper, we only detect one-dimensional Jordan blocks. 

We now ask what happens for different open system effects corresponding to one-dimensional Jordan blocks. We 
focus on the case where we have s^g/ = 0, J being an index set, and Re(si^/) < 0. In the t ^ oo limit, the matrix 
e"'* tends to the matrix P' , which projects onto the subspace corresponding to /. Thus, in this Hmit, e^"^* = Se'^*S~^ 
tends to the projector SP'S~^ and the asymptotic coherence vector a;'°°' depends on the initial coherence vector xo 
according to 

a;(°°) = Pxa - Ar^b. (35) 

On the other hand, if / is the empty set, i.e., Re(si) < Vi, we have e^^* — Se'^^S^^ — > when t ^ go. In this case, 
the asymptotic state is determined by 

a;(°°) = -M-^b, (36) 

thus being independent of Xq. In this paper / is the empty set, we therefore proceed to analyze the latter asymptotic 
state in the next section. 



C. Susceptibility for the asymptotic state 

Let us now examine the susceptibility for a small probe field that interacts with a medium consisting of a collection 
of A-systems that in turn are affected by the open system channels considered in Sec. IIII Al The susceptibility of a 
medium is related to the asymptotic state of the atomic coherence p^^^ = xf^'^ + ix'c^'^ , which can be calculated by 
using Eqs. (021), and see., e.g., Ref. Explicitly, 

X = ^pir'M, (37) 

with K = iV|(a|/i|6)p/(2eo?i), and we recall that the Rabi frequency is defined in a non-standard way with an extra 
factor of 2 for notational convenience: Vli = {a\fi\i) ■ E/(2fi,), i = b,c. The susceptibility defined in Eq. l(37j) is correct 
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0.0 0.1 0.2 0.3 0.4 0.5 0.0 0.1 0.2 0.3 0.4 0.5 



l"c|/7 IHcI/t 

FIG. 2: (Color online) In (a) and (b) we show the real and imaginary part of the susceptibility. In (c) and (d) the slow-down 
factor (group index) Ug and absorption a at a given wave length A, on two-photon resonance are shown as a function of the 
control field \Q,c\. The thin black line shows the case of no open system interaction on Vbc- The thick gray line shows the 
dephasing channel and the thick black line shows the depolarization channel. The remaining lines describe three versions of 
the amplitude damping channel, the wide-dashed (blue) line shows fe <— c damping channel, the dashed (red) line shows c <— fe 
damping channel and the dash-dotted (green) line shows the incoherent sum of the two damping channels with equal rate. In 
all plots the proportionality factor k defined below Eq. I|37ll is scaled to 1. In (a) and (b) the numerical values are \Q.c\ ~ 0.16, 
7 = 1, A = 0, and all non-zero channel rates are set equal to 0.1. In (c) and (d) the channel rates are normalized so that the 
slow-down factors are matched. 



in the limit of a small probe field, since in this case pab is linearly dependent on fJ;,. Furthermore, we recall that the 
real part of the susceptibility is related to the phase evolution of the optical field, while the imaginary part describes 
the absorption or gain of the field. 

The open system channels defined by the Lindblad vectors in Eqs. (|22l) . l(25l) . and l(27|) . change the properties of 
the A-system in a nontrivial way. The susceptibility in the general case can be calculated from Eq. ((36l) and we find 





2r*Re(7c)(7 + 277J + 77:(2r?_ - 


+ 7(2r7_ +?/jr*|7,|2 


{\n,\^ + T*X) 


|r!,|2Re(7,)(27 - 2r^_ -I- 3r; J -I- 277?^ |7,|2 





where the parameters are defined in Eq. (|3T|) . This is a general expression for the susceptibility when the medium 
interacts with a Markovian environment through the open system channels. A general observation is that this 
expression does not depend on the optical phases even though the dark state does. 

Below, the susceptibility due to the different open system channels are studied one by one. Analytical expressions are 
retrieved from Eq. l(38|) as special cases by setting all channel rates except one to zero. To facilitate easier evaluation of 
the characteristics of the susceptibility due to the different open system channels we plot the corresponding quantities 
in Figs.[2l[a) andEUb) using the same channel rate for all open system channels. Since the susceptibility is a complex 
function, the real and imaginary parts, are plotted separately. 
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If we put all channel rates to zero except rjz we find the susceptibility due to dephasing. Accordingly 

(S + iVz) 



Xdephase 



(39) 



This expression for the susceptibility is frequently encountered in the literature, e.g., Ref. [20|. 

If we instead put 2?7_ = ?7+ = 776c 7^ we find the susceptibility of a medium that effectively changes the population 
within the subspace Vbc and is symbolized accordingly: b ^ c. We have, 



(S + i?76c/4) 



'(5 + jW4)('5 + A + i7)- |0elV2' 



(40) 



This expression for the susceptibility is also frequently encountered in the literature, e.g., [27|. It is interesting to note 
the similarity of this expression with Eq. (|39|) . although the open system channels are quite different. Nevertheless, 
in Figs.[2ta) andElJb) it can be seen that the phase response of Xb^c is larger than for Xdophasc- This means that a 
medium affected by pure dephasing cannot slow down light as efficiently as compared to a medium affected by the 
b ^ c channel, assuming that the channel rates are equal. 

We now consider the amplitude damping channel in the other direction, c <— 6. The analytical expression of the 
susceptibility is given by 



Xc^b 



2|7(4z(5 + rjcb) + r]cb 


77cf,-2i((5 + A) 






{{r]cb--iiS) 


2^7 + 2{5 + A) + irjcb 


+ 4i|f^cP}{ 


72 + ((5 - 


- A)2 


Vcb + il + 2Vcb)\^c\^} 



(41) 



This channel shows a different behavior compared to all other considered channels, it shows gain where the other 
channels show absorption. That is, a slowly propagating light pulse will be amplified provided that the open system 
channel c ^ b can be established. 

Let us now consider the channel that is defined by rjpe = ijbc = Vcb, i-e., the amplitude damping channels in both 
directions. This open system channel may be important in experimental situations where, e.g., atoms collide and the 
populations of the states undergo sudden incoherent exchange. We denote this as the population exchange channel, 
'popex' for short, and the susceptibility Xpopox- Hence, 



-777pe(2(5 + iripe) 


'4{S - A)2 + (27 + j/pe)"' 


+ 47 


- 27(5 + (-2,5 + A)77pe 




\^{2S + iT]pe) 


2^7 + 2(5 + A) + ir]pe 


-2|r!cP}{7?ype 


A{S 


- A)2 


+ (27 + Vpef 


+ (27 + 


ripe){l + 3?7pe) ficpj 



Xpo 

<J (2,5 + ir^pe) 2^7 + 2(5 + A) + ii^p^ - 2|17c|2 f- 1 7rype 4(5 - A)^ + (27 + r/pe)^ + (27 + 77pe)(7 + Sr/pe) 

(42)' 

As expected, the absorption profile is located in between the profiles of Xb^c and Xc*-6, as can be seen in Fig. [2l[b). 
To first order, the absorption profiles of the different amplitude damping channels is shifted by the parameter 77 — 
(j]bc — Vch)/2. However, the phase response, i.e., Re(x6^c) is not affected by this shift, and it is approximately the 
same for all considered amplitude damping channels. 

The susceptibility associated with the isotropic depolarization channel (77 ^ 0) is given by 



Xdcpol — 1^ 



2i-f^r]{3d + 2iT]) 


4(5 - A)2 + (27 + 77)2 


+ 3 


7(65 + 2ir]) + (55 - 2A + 777)77 




{(35 + 2i77) 


27-2i(5 + A) + ?/ 


+ 3i|17,|2||2777 


[4(5 - A)2 + (27 + r;)2 


+ 3(27 + 77)(7 + 277)|r!,|2} 



(43) 



In experiments, e.g., where one uses laser beams and hot atomic gases, this open system interaction should be relevant 
as it corresponds to coherently prepared atoms leaving the laser beam, replaced with atoms in a completely mixed 
quantum state. The susceptibility is plotted in Fig. ^ and it can be seen that the phase response, Re(xdcpoi), is 
between the corresponding curves for the dephasing and the amplitude damping channels. Furthermore, the absorption 
curve, Im(xdcpoi), is close to the absorption profile of the depolarization channel, but significantly smaller than the 
corresponding curves for the dephasing and b ^ c channels. With respect to absorption and phase response, the 
depolarization channel is therefore not as harmful as the dephasing channel. 

When there is no open system interaction on Vbc we retrieve the ideal susceptibility. 



25 



Xidcal = -K 



25(5 + A + n)- 



(44) 



10 



included in Fig. ([2]) as a reference. 

Often one is interested in achieving the slowest possible group velocity in a medium consisting of A-systems. Thus, 
we calculate and present the slow-down factor as a function of the control field, this is plotted in Fig.[2Kc). The slow- 
down factor, or group index, is proportional to the derivative of the real part of the susceptibility, Ug = LudgKe (x) /2 at 
angular frequency ui. In Fig.[2]Jd) we present the absorption on two-photon resonance as a function of the control field, 
where the absorption coefficient is defined as a = (27r/A)Im To compare the absorption for a given slow-down 
factor at (5 = A = we normaHze the channel rates. Approximating dsXk for small channel rates in Vbc gives 



dsRe (Xdophasing) 

dsRe ixb^c) 
dsRe {xc<-h) 
dsRe (xpopcx) 
dsRe (xdopoi) 



2k 



S=A=0 



5=A=0 



2k 
2k 
2k 



S=A=0 



<5=A=0 



2k 



(45) 



thus, the normaHzed rates are r/z = ?76c/4 = r]cfc/4 = ?7pe/4 = 2ri/3. In Fig.jSJd) we see that the different open system 
channels give different absorption profiles, even though the slow-down factors are equal in this regime. Furthermore, 
it is interesting to note that the dephasing and the b ^ c amplitude damping channels, with normalized channel rates, 
have identical absorption characteristics. We see that the isotropic depolarization and population exchange channels 
are quite similar in Fig. [2][d), but more essential, they have a lower absorption than the dephasing and the b ^ c 
amplitude damping channels. As already mentioned, the damping c ^ 6 is special in that it yields a gain where the 
other channels yield absorption of the probe field. 



IV. CONCLUSIONS 



We analyze a three- level A-system interacting with a Markovian environment. The Lindblad master equation is 
reformulated using the Gell-Mann matrices as a basis. We show that we can rewrite this in a convenient vector 
formalism that is independent of the specific matrix representation of the SU(3) generators. An evolution matrix 
corresponding to this vector form is worked out in detail and the dynamics is written as an inhomogeneous matrix 
differential equation. The general solution of this equation can be given in terms of the Jordan normal form. The 
Jordan blocks of the evolution matrix are one-dimensional and the real part of the spectrum is negative. From this 
follows that the asymptotic solution in all cases we are interested in are independent of the initial states and given by 
the inverse of the evolution matrix multiplied by a vector that is given by the open system interaction. 

The Hnear optical response of a quantum system is given by the asymptotic steady state solution of the evolution 
equation. Assuming a weak probe field and a strong control field, i.e., the familiar settings of electromagnetically 
induced transparency (EIT), we find an explicit and closed form for the susceptibility for general open system dy- 
namics, including pure dephasing, depolarization, and three different forms of amplitude damping channels acting on 
the subspace of the ground states of the A-system. 

The open system effects introduce a non-trivial behavior of the linear optical response of the A-system. For 
example, we show that the isotropic depolarization channel is less detrimental to the induced transparency and the 
corresponding slow-down effect, than is pure dephasing. Furthermore, an amplitude damping channel can lead to 
reduced absorption and even gain for the probe field. 

It is important to know how and why the coherence vector is changing in direction and length, not only in areas 
of quantum information and computation, but also in high precision spectroscopy, metrology and photonics. The 
results of this paper shows the usefulness of the coherence vector formaHsm for experimental and theoretical studies 
of light-matter coherence phenomena. 
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